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1  Summary  of  Research  Accomplishments 

The  majority  of  the  research  focussed  on  developing  new  anti-jam  and  code  synchronization  algorithms 
for  GPS  receivers  (funding  from  NAWCWPNS).  In  addition,  support  for  a  graduate  student  (Craig 
Agate)  was  provided  by  BMDO  for  part  of  the  period  9/30/98  to  12/30/98,  during  which  time  he 
started  writing  his  dissertation.  We  anticipate  that  the  final  dissertation,  [1],  describing  work  completed 
under  both  this  contract,  and  the  previous  ONR  contract  N00014-95-1-1213,  will  be  completed  in  the 
Summer  of  2000.  The  following  specific  algorithms  were  developed  in  this  reporting  period. 

1.  Modulated  Lapped  Transform  (MLT)  for  GPS  C/A  code.  The  MLT  is  related  to  Fourier  and 
wavelet  transforms,  and  generates  a  time-frequency  decomposition  of  the  received  signal. 

2.  Extended  Kalman  Filter  (EKF)  for  code  tracking  using  the  MLT  (MLT-EKF).  This  algorithm 
jointly  estimates  the  code  delay,  jammer,  and  received  code  amplitude.  The  jammer  estimate  is 
subtracted  from  the  received  signal  as  part  of  the  EKF  algorithm. 

3.  Generalized  Likelihood  Ratio  Test  (GLRT)  for  C/A  code  acquisition.  The  EKF  requires  prior 
knowledge  of  the  code  delay  to  within  half  a  chip.  To  provide  initial  acquisition,  a  GLRT-based 
algorithm  was  developed  which  provides  constant  false  alarm  rate  (CFAR.) 

The  GPS  code  tracking/anti-jam  algorithms  based  on  modulated  lapped  transforms  (MLT)  are  now 
summarized. 


2  Transform-Domain  Model  for  the  GPS  Signal/Interferer 

In  [2]  we  developed  a  vector  model  for  the  received  GPS  C/A  code  signal  using  the  MLT.  The  discrete¬ 
time  received  signal  is  given  by 

r{k)  =  a(kTs)PNip(kTs  -  r{kTs))  cos(2ttFcA:Ts  +  6)  +  j(k)  +  n(Jfc),  (1) 

where  Ts  is  the  intermediate  frequency  (IF)  sampling  interval.  The  jammer  j(k)  and  additive  thermal 
noise  are  real- valued,  bandpass  processes,  as  is  the  modulated  C/A  code  PN[p(t ).  Note  that  the  signal 
must  be  real-valued  to  use  the  MLT.  The  sequence  PNip(t)  is  a  filtered  pseudo-random  binary  waveform, 
with  chip  rate  1.032  mcps. 
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The  overall  estimation  problem  is  to  (a)  acquire  and  track  the  code  delay  r(kTs),  while  simultane¬ 
ously  (b)  rejecting  the  jammer  j(k). 

A  general  form  of  the  jammer  was  adopted  in  [2],  encompassing  both  sinusoidal  and  chirp  waveforms. 
Thus, 

j{k)  =  \fj cos(27t((Fc  +  Fj)kTs  +  Fj(kTa  -  mTswp)2)),  (2) 

where  Fj  is  the  jammer  center  frequency,  and  2 Fj  is  the  sweep  rate  in  Hz. /sec. 

In  order  to  apply  the  MLT,  the  received  samples  r{k)  are  grouped  into  vectors  of  length  4 Nc  samples, 
where  Nc  is  the  number  of  chips  per  vector,  and  1  /Ts  =  4 /Tc  is  the  sampling  rate.  The  received  vectors 
are  then  compactly  written  as 

r (k)  =  S(k,  T(k))a{k)  +  j(fc)  +  n(k),  (3) 

where 

r(Ar)  =  [r(kM),r(kM  +  1), . . .  ,r({k  +  1  )M  -  1)]T  (4) 

The  signal  matrix  S (k,r{k))  €  TlMx2  is  defined  in  [2],  Note  that  the  received  amplitude  has  been 
separated  into  its  inphase  and  quadrature  components,  yielding 

a (k)  =  [ a(kMTs )  cos (9),a(kMTs)  sin(0)]T  .  (5) 

The  MLT  is  defined  by  a  transformation  matrix  T  €  7lMx2M  [2],  The  transformed  received  vector 
used  to  develop  the  code  acquisition  and  tracking  algorithms  is  finally  defined  by 

r(k)  =  T[r(fc  -  1)T, r(fc)T]T.  (6) 

Note  that  T  transforms  a  2 M  length  vector  to  an  M  length  one,  and  furthermore,  is  overlapped,  in  that 
both  r(jfc)  and  r(k  +  1)  depend  on  the  same  time-domain  vector  r(k).  The  transform-domain  received 
signal  model  is  then 

r(k)  =  S(k ,  r(k))a(k)  +  n'(fc),  (7) 

where  n'(fc)  denotes  the  transformed  noise  plus  jammer. 

The  major  advantage  of  the  lapped  transform  is  that  the  correlated  time-domain  jammer  plus 
noise,  j{k)  +  n(k )  is  now  converted  to  an  uncorrelated  vector  sequence  n'{k).  This  greatly  simplifies 
development  of  the  EKF  algorithm,  as  described  in  the  next  section. 


3  Modulated  Lapped  Transform-domain  EKF  and  GLRT-based 
Acquisition 

In  order  to  apply  the  EKF,  a  state  vector  comprising  the  delay  and  amplitude  is  first  defined  by 

x(fc)  =  [r(/c),a(fc)T]T  (8) 

An  EKF  form  of  the  measurement  equation  is  then  given  by  [2] 

r(k)  =  h(k,x(k))  +n'{k),  (9) 

where  the  nonlinear  measurement  function  is 

h(fc,x(fc))  =  S(fc,r(fc))a(fc).  •  (10) 

The  EKF  can  be  directly  applied  to  the  measurement  model  (9),  and  the  resulting  algorithm  is 
simplified  since  n'{k)  is  approximately  white  Gaussian  noise.  However,  the  covariance  matrix 

R(lfc)  =  E{n/{k)n/{k)T} 

is  unknown  due  to  the  presence  of  the  jammer.  A  method  based  on  maximum-likelihood  is  presented 
in  [2]  for  estimating  this  matrix. 
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Initialize  ACQJFLAG  =  0 
For  l  =  0, 1, . . . 

Input  next  received  vector  r (l). 

Compute  MLT  vector. 

L(l)  =  T[r (/  -  l)Tr(/)T]T 
Update  noise  covariance  matrix  estimates 

*<>(0  =  tz.\i+iL(1)l(1)t  +  R0(7  -  1) 

Ri(0  =  T3prr[r(0  -  S(Z,f(Z|Z  -  l)a(Z|Z  -  1  )][r(Z)  -  S(l, t{1\1 —  l)a(Z|Z  -  1)]T 

+At^tRo(/-1) 

Update  error  covariance  matrix. 

pm-1  =  p(i\i  - 1)-1 + h(otRi  (o_ih(/) 

Compute  EKF  measurement  update. 

x(/|/)  =  x(Z|Z  -  1)  +  P(Z|Z)H(Z)TR1(Z)_1[r(Z)  -  S(Z,f(Z|Z  -  l)a(Z|Z  -  1)] 
Compute  one-step  predictions. 
x(Z  +  1|Z)  =  Fx(Z|Z) 

P(Z  +  1|Z)  =  FP(Z|Z)FT  +  Q 
Test  for  acquisition  using  GLRT. 

If  |Ro(Z)|/|R1  (Z)J  >  Th  then 
ACQ.FLAG  =  1 
r(Z  +  l|/)  =  r(/  +  l|0  +  ^BC 
Else 

ACQ_FLAG  =  0 

Next  l 


Table  1:  MLT-EKF  Algorithm  Summary. 


The  acquisition  algorithm  is  based  on  the  GLRT.  Let  hypothesis  H0  and  Hi  correspond  to  not 
acquired  and  acquired,  respectively.  Then  the  measurement  models  are 


H0  :  r(k)  =  ^(k)  (n) 

Hi  :  r(k)  =  S(k,r(k\k  -  l)a(fc|A:  -  1)  -l-n^A:). 


That  is,  if  acquisition  is  occurring,  the  received  vector  consists  of  the  estimated  signal  vector,  thermal 
noise,  and  jammer.  However,  if  the  estimate  f(k\k  -  1)  is  far  from  r(Jfe),  it  is  assumed  that  the  C/A 
signal  component  is  effectively  a  white  noise.  Hence,  under  Ho,  the  received  signal  is  a  Gaussian  noise 
vector  n^ik). 

The  GLRT  is  defined  by 


* 

maxRnitoPfcffll/^R)  > 

maxRnf=o  P(H(Z)  Iffo.R)  < 

#0 


(12) 


where  Th  is  a  threshold,  and  R  corresponds  to  the  unknown  covariance  matrices  of  ^(k)  or  nx{k). 

Acquisition  is  performed  as  follows:  If  the  GLRT  is  greater  than  the  threshold,  acquisition  is  declared, 
and  the  MLT-EKF  updates  normally.  Otherwise,  the  delay  estimate  is  updated  by  a  fixed  valued  rinc’ 
to  step  through  the  code  delay  uncertainty. 

The  final  MLT-EKF  algorithm,  with  GLRT-based  acquisition  is  summarized  in  the  following  table 
from  [2].  6 


4  Example  Results 

An  example  result  from  [2]  is  shown  in  Figures  1  and  2.  In  Fig.  1,  a  swept  tone  jammer,  with  2 Fj  =  2000 
mHz. /sec.  interferes  with  the  C/A  code  waveform,  received  with  an  SNR  of  -5  dB.  Again,  without  AJ, 
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Delay  Estimate  (Chips) 
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Figure  1:  C/A  code  acquisition  and  tracking,  with  and  without  interference  rejection  for  a  swept  tone 
jammer,  2000  mHz.  sweep  rate  . 


the  MLT-EKF  algorithm  fails  to  acquire.  The  MLT  bin  powers  (diagonal  elements  of  R (k))  are  shown 
in  Figure  2  for  this  swept  tone  jammer. 


5  Conclusions 

New  algorithms  for  jammer  suppression  and  joint  C/A  code  acquisition/tracking  were  developed  for 
GPS  receivers.  The  modulated  lapped  transform  was  employed  to  decorrelate  the  jammer,  and  provide 
better  isolation  of  time/frequency  content.  As  a  result,  the  overall  algorithm  was  able  to  reject  both 
tone  and  chirp  jammers  with  high  sweep  rates. 
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Figure  2:  Averaged  output  power  of  MLT  bins  for  the  swept  tone  jammer. 
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